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INTRODUCTION 


The purpose of this vork is to review several methods "by which a 
magnetic field B in space can he represented, with particular 
attention to problems of the observed geomagnetic field.Time depen- 
dence will he assumed to be negligible and five main classes of repre- 
sentation will be described, as follows: 


(1) Representation by a vector potential A : 

B = V A (l) 

This is a general form to which other forms noted here can 
always be converted* 

(2) • Representation by a scalar potential 7 : 

B = - V 7 (2) 

This representation is available- only in current-free regions and 
it is particularly appropriate to the main geomagnetic field, where 
the expansion of 7 in spherical harmonics is of considerable 
interest* 

(5) Representation by orthogonal vectors, in particular those 

related to spherical coordinates 


This representation is related to spherical vector harmonicsj 
it has been used in dynamo theory and in a variety of problems. 
As will be seen, it is also useful in devising models of the 
earth's magnetosphere field. 
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(4) Representation 'by Euler potentials 

B = Vot X Vp (4) 

This is the only representation vhich includes explicit infor- 
mation about the configuration of magnetic field lines, hut its 
nonlinear character m^es its derivation difficult. In the 
earth's magnetic field c>d is related to Mcllwain's L parameter 
■which is useful in the study of the motion of trapped particles. 

(5) Local representations, in which B is expanded around its 
value at some reference point 

B= BQ+r*VBQ+ |rr:VVBQ+... (5) 

Such expansions find use in the theory of guiding center 
motion and in describing the vicinity of neutral points at 
which B vanishes. 


The present discussion is not intended to be self-contained: matters 
on which recent reviews exist in the general literature will only be 
briefly described with references directing the leader to more elaborate 
treatments. References will also be given to articles in which various 
mathematical tools described here are employed or mentioned, but the 
mentioning of such articles does not always imply concurrence with all 
conclusions expressed there, nor is the review of the literature meant 
to be complete. In some cases review papers are cited mainly in order 
that they may provide the reader with a guide to earlier work not 
mentioned here. 
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(1) THE VECTOR POTEHTIAL A 

The representation (l) is the most general one and others may be 
reduced to it. For instance, (3) leads naturally to 

^ E ( 6 ) 

■while (if) gives 

A = (7) 

■which is orthogonal to B . The representation (2) also has an equivalent 
form (l) and if y is expanded in spherical harmonics an equivalent 
representation of the form ( 3 ) is readily obtained (as ■will be shown), 
allowing (6) to be used* All these choices of A are indeterminate 
within the gradient of an arbitrary scalar 'Vq , since the addition 
of such a gradient to A does not affect (l) . 

The representation (l) is part of the standard treatment found in 
practically all texts on classical electrodynamics and will therefore 
not be discussed in detail* Its main usefulness, in that form, is that if 
the current density j creating the field is everywhere known, A is 
readily" derived (assuming B to be a vacuum field) by volume Integration 

j(r' ) 

A(r) = (hqAJT) J dV' (8) 

In engineering applications j is often given by the circuit 
geometry and A is readily calculated (often, B is directly derived in 
such cases from the Biert-Savart formula). In space, on the other hand, 

J is generally a derived quantity obtained only after B is given, 
often with low accuracy, so that (8) is used relatively infrequently. 

One important case In space science ■where j given occurs in 
regions where it is inferred from the observed fluxes of charged particles 
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Northro p, 1963 , eq. 4-1 ; Lon^lre , 1963 ] . One such region 
occurs in the vicinity of the earthy where trapped particles support 
a current density often called (not entirely accurately) the n'ri g current - 
The magnetic field in these regions is generally imown from independent 
observations and its comparison to j forms an interesting check, often 
enabling one to deduce which part of an observed particle population 
contributes most to the magnetic perturbation- Many analyses of the 
ring current obtain B directly from j via the Biot-Savart formula 
Akasofu and Chapman , 1961 ; Schield , I 969 J = However, general 
formulas for A , derived by (8) for configurations used in calculations 
of the effects of a ring current, have been published by Kendall et al , 

[ 1266] . 

Another case in which distributions of current density are used 
as the source of magnetic fields in space occurs when such distributions 
are introduced as the input data of theoretical models- The configuration 
of the current in such cases may be selected on physical grounds - e-g. 
field-aligned currents Bonnevier et al. , I 97 Q ; Crooker and Siscoe, 

] or the geomagnetic tail sheet current [~ W.i Hiams and Mead , I 963 - 

or else it (or part of it) may be represented by a generally expanded 
function of position. The advantage of using J In the latter case is 
that this assures the vanishing of 7 *B , which may have been the 
reasoning of Olson 1974 ; Olson and Pfitzer , 1974 ] . Later on in 
this review, however, it will be shown that similar properties may be 
obtained by 'simpler means. In most such applications, the Biot-Savart 
formula is used and B is derived directly. 



- 6 - 


(2) THE SCALAR PCOEriTIAL 7 


(a) IT5TRODUCTION 


In current-free regions 7 X B vanishes and it is possible to 
represent B by a scalar potentisil 

B = - V7 (9) 

Because = 0 , 7 is harmonic 

= 0 (10) 

and may be expanded in a variety of vays appropriate for harmonic functions. 
The properties of harmonic functions and of y in particular are dis- 
cussed in most texts of classical electrodynamics and for this reason 
the discussion here vill be confined to applications relevant to the 
geomagnetic field. 


(b ) SPHERICAL HARMONICS 


The "main" geomagnetic field - i.e. that part of the field observed 
at or above the earth’s surface which is caused by currents in the 
earth’s interior - accounts for about 99 °/o of the field observed 
at ground level and is readily expanded in external spherical 
harmonic functions 

7 = a ^ (a/r)^^^ 1^(0) (g™ cos mi* + sin mi ) (ll) 

_ n n n 

n=l m=0 

where a is the earth's radius and are associated Legendre functions. 

Attention should be given to the factor a outside the summation, inten- 
ded to give the coefficients and h^ the dimensions of magnetic 
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field : as a consequence such coefficients are generally given in units 
of gausses or gammas ( one gamma =17= 10"^ gauss ), The form given in 
(ll)j involving real coefficients, is preferred in geomagnetic research 
over the one used in mathematical treatments, vhere the last factor in 
each term of (ll) is expressed in terms of exp (imM’) , with the 

summation over m extending from - n to n and the coefficients 
} "being complex conjugates ( a related representation 
uses spherical harmonic functions proportional to P^(0) exp im'f ) , 

Properties and details of this expansion are reviewed in many texts 
and articles [e.g., Chapman and Bartels , 1940 ; Heppner , I965 ; Kaula , 

1968 ; Stacey , 1969 ; 21muda , 1973 ] 

In using the spherical harmonic expansion of 7 note should he taken 
of the choice of normalization, for several definitions of Legendre 
functions, differing by constant factors, are currently employed. 


Three main choices of this kind are described by Chapman and Bartels 
£ 1940 ] ; in the notation used by them (adopted from Adolf Schmidt) 
Legendre functions are denoted P for mathematical normalization, 

n, m rn 

for gaussian normalization and for Schmidt normalization. The 

same notation will be used here but it is by no means a general one 
and frequently authors denote Legendre functions by regardless of 

the normalization which they are using. 


In the above notation, if p = cos0 

^n 0^*^^ " ^ (-l)''/2'' nl ) d'^/dp^ (1 - p^)"" 

(Legendre polynomials) 

= (1 - dXy ap” 

= [(n - m)J / 1.5. .. (an - l)] 

= [(n - m): nl a“ / (2n);] P (n) 

t n,m 


( 12 ) 

(13) 


(U) 
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p 

n 

^n,0 


(15) 

II 

[ 2 (n - 

m) V(n + m)J ] P 

J n,m 

B 

0 

p° = 
n 

[1.3.5. 

... *(2n - 1 )/ nl'l 

(16) 

II 


... *(2n - 1)/[| (n t m) 

i (n - m)i] 




(m ^ 0) 


The coefficients of (ll) generally have their indices arrayed in 
the same manner as the Legendre functions to which they are attached# In 
transforming (ll) from one normalization to another, whenever a Legendre 
function is multiplied by seme factor^ the corresponding harmonic coeffi- 
cient has to he divided hy the same factor, ensuring that y does not 
vary# 

In the literature coefficients of the e^^ansion of y are generally 
listed for either Schmidt or gaussian normalization: for easy identifi- 
cation the convention exists to reverse all coefficient signs for gaussian 
normalization, so that the axial dipole coefficient is negative for 

Schmidt-normalized potentials hut positive for gaussian ones. In Schmidt 
normalization the magnitude of terms is roughly of the order of their 
contribution to the field: their gradual decrease with growing n reflects 
the relative preponderance (near the earth's surface) of the contribution 
of low-order harmonics, and for any n they do not vary systematically 
with m , 

In gaussian normalization the magnitude of terms decreases more 

graduall/ with n and for any n the terns having low m tend to be 

larger. In computer applications the derivation of P is, however, 

n • ni 
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much simpler than that of , because it can be handled by means of a 
recursion relation free from any irrational factors, connecting terms 
vith the same m f Cain et al . , I967 ] : 


P (0) = sin^9 

m,m 


P (0) = cos 9 sin°^ 0 

m+l,ra 


(17) 


P T (0) 
nfl,m 


cos 0 P 


n,m 


(n -t* m) (n - m) 
(2n *t- l}(2n - 1) 


■n-l,m 


The derivatives required for the calculation of are best found 
by a recursion relation based on the derivatives of the above relations, 
starting -with 


dP /d0 = 
m.m 


. ra-l« 
m COS0 sin 9 


(18) 


Normalizations other than those described also exist: a comparative 
list of 8 choices from I8 choices has been compiled by Kaula 

r 196^ ; Table ij . 


(c) GEOMAGNETIC MOKELS 


Expansions of the form (ll) have been used for the geomagnetic field 
since the time of Gauss 1839 ] • Typically, model expansions use 
harmonics vith n ranging up to n 9 with n between 6 and 10 
Since the number of independent coefficients in such an expansion is 
(n + 1)^-1 , one often finds authors speaking about models with 

48, 63 , 80 , 99 or 120 coefficients. 
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A very comprehensive review of early analyses of the main field, has 
been compiled hy McDonald and Gunst ^ I967 ] and some early results are 
also given hy Chapman and Bartels ig^oj - More recent field models 
have been listed by Heppner ^ ^ 9 ^ 3 ] Cgaussian normalization) and by 
Ve Stine 196 o j ( Schmidt-normalized) . Cain 197 l ] reviewed some other 

recent work and analyzed problems which arise in connection with 
field mapping by satellite. An International Geomagnetic Reference Field 
IGRF 1965.0 was derived and published by lAGA Commission 2 , Working 
Group k- ^ 1969 ] and contains 80 terms. 

Most of these models take into account the slow "secular" variation 
of the field by assuming a linear dependence of the coefficients h^) 

of the form 

a = At (19) 

where a^ is the value of the coefficient at some initial time ( e.g. 

the beginning of I965 for the IGRF model mentioned before) and t is 

in years. Some models - for example that of Cain et al , [ 1967 ] - also 

2 

include correction terms proportional to t • 

(d) CURRENT-FREE MODELS OF THE MAGNE?TOSPHERB 


The expansion (ll) consists solely of "internal harmonics" representing 
a field which originates inside the earth and vanishes at infinity. If 
sources external to the earth were' also included their contribution to 7 
could be represented by a series of "external harmonics" with positive 
powers of r 

^ m=n 

7 ' = a ^ ^ ^ (r/a)^ t ^ ^n ^ J 

n=l m=0 


Models of this form have been proposed for the magnetospheric field 
by Mead [12a] and by Midgeley [12^] , In both cases the earth's 
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field "was approximated "by the axial dipole given "by in (ll) j "which 

is orthogonal to the earth-sun line_, and the frame of coordinates for 
both the dipole component and the expansion (20) is fixed with respect 
to the sun and the plane of the ecliptic* 

In more complicated models internal and external fields are usually 

expressed in different frames of reference, since the internal field alone 

co-rotates with the earth. If one neglects all internal harmonics except 

for the three dipole coefficients gj and hj , the model will 

depend (when external conditions are fixed) only on one parameter which 

varies in time - the angle between the earth's dipole axis and the 

0 o 

eart-sun line, which varies by about +55 around 9 ^ , depending on 

the season of the year and the rotation of the earth* For any given 
value of this angle such a model is symmetrical about the plane containing 
the dipole axis and the earth-sun line and it is conveniently described 
in solar geomagnetic coordinates F Olson , 1969 (where the caption of Fig* 1 
should be ignored)j Bussell , 197l > Burch and Janetzke , 1974 ^ with the 
origin at the earth's center, the x-axis pointing sunward, the x-z plane 
containing the dipole axis and the z-axis pointing into the northern 
hemisphere • 

Current-free models of this kind are not meant to represent exter- 
nal fields observed on the earth's surface - indeed, the surface effects 
predicted by such models are far too small to account for the observed 
daily variation [ Mead , 1964 ] * Instead, they are intended to approximate 
the large-scale field of the outer magnetosphere. In a qualitative way 
they seem to agree with observations - in particular, they display a 
sharp boundary on the sunward side, limiting the region of field lines 
connected to the earth and corresponding to the observed magnetopause. 

The most significant non-dipole coefficients seem to be and gi 

and further improvement has been obtained \ Williams and Mead , 1965 ] 
by adding a current sheet across the median plane of the geomagnetic 


tail. 
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Another simple current-free model approximating the external magneto- 
sphere is the image dipole models here the effects of the external field 
are approximated hy a dipole parallel to the z-axis hut located some 
distance sunward from the earth Hones , 1963 ; Taylor and Hones , I 965 ; 
Forhes and Speiser , 1971 ^ » This representation, which was inspired hy 

the theoi^” of the magnetic storm developed hy Chapman and Ferraro ^ I93lj , 
has two adjustable parameters - the distance to the "image dipole" and 
its magnetic moment, which generally exceeds that of the earth by a 
considerable factor.. For instance, in the work of Taylor and Hones I 963 J 
the image dipole is 28 times stronger than the axial dipole of the 
earth (other internal terms are ignored) and is placed earth radii 
sunward of the earth, at a point which is outside the "magnetopause" 
and thus beyond the region in which the model is valid; the model also 
includes a sheet current in the geomagnetic tail region, somewhat 
similar to the sheet current introduced by Williams and Mead ^ 1963 J o 
In all models with image dipoles the scalar and vector potentials are 
readily found by superimposing the contribution of the two dipoles, and 
it is usually best to leave them in this form and not expand 7 in 
spherical harmonics. 


(e) CHANGE OF COORDINATES 


In general the expansion (ll) is given in spherical coordinates with 
the origin at the earth's center and the z-axis aligned with the terres- 
trial rotation axis. For some applications it is useful to rotate the 
z axis so that it coincides with the dipole axis, giving the so-called 
"dipole coordinates" ^ Chapman , 1963] « If 7 is expanded in dipole 

coordinates both gj and h^ vanish and the entire dipole contribution 
is contained in the axial dipole term , 

The use of dipole coordinates is especially advisable when it is 
desired to enhance the pieponderance of the axial dipole component over 
other harmonic terms. For instance, geomagnetic Euler potentials 
(described later) may be obtained by a perturbation calculation in which 
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the zero-order term is the contribution of the dipole field; it is then 
clearly advantageous to use dipole coordinates vhich make this term 
relatively large and reduce the remaining "perturbation". 

Mathematical formulas exist for transforming (ll) from geographic to 
dipole coordinates and they have the virtue that harmonics of a given 
lower index n contribute in the new coordinates only to harmonic terms 
■with the same n , Thus each group of 2n tl coefficients corresponding 
to a given n transforms independently and (say) a-model with 99 coeffi 
dents ( n ^ 9 ) will be represented with full accuracy by the same 
number of terms in the new coordinates (to be precise^ there will only 
be 97 independent terms in the new expansion, since and h^ vanish) 

formulas for such a transformation have been discussed by [1969 j 

and are concisely given by Slater [1260] . 

Where high-speed computers are available it is often simpler to 
apply a %rute force" transformation as follows. Suppose a 99-term 
expansion of y Is given: one then begins by deriving its value 
from the given expansion (in geographical coordinates) at 97 points 
scattered over the surface of the earth. Next one transforms the 
coordinates of each of the selected points to (spherical) dipole 
coordinates in which the z axis is antiparallel to the vector having 
cartesian components ( g^, ) • In the new coordinates the expan- 

sion of y at each of the given points is derived, with the new 
expansion coefficients (g“', h“’ ) entering as 97 unknown quantities 
(t-wo of them vanish and are not counted). This produces 97 equations in 
97 unknowns, which are now solved to give the new coefficients, A repeti- 
tion based on another set of points is recommended as a check on accuracy 
and as a precaution against ill-conditioned sets of equations, but the 
resulting set of coefficients is usually accurate enough for normal use. 
An example of this method is given by Stern [1971 ] who included in the 
transformation not only the time derivatives of (I9) but also the slow 
secular variation of the new coordinate axes,- 
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If the origin of the coordinates is also alloved to vary "we obtain 
eccentric dipole (or ’’offset dipole") coordinates* The choice of an 
arbitrary origin introduces 3 more adjustable parameters and can be 
used, for instance j to assure the vanishing not only of and hj 

but of gg j gg and as well |^ Bartels , 1936 ; Chapman and Bartels , 

19lj.O j . The z axis of eccentric dipole coordinates should parallel 
that of dipole coordinates to assure the vanishing of the two off-axis 
terms with n = 1 


Unfortunately, a given expansion in geographic or dipole coordinates 
(e.g. with 99 coefficients) is no longer accurately represented by the 
same number of coefficients (or in general, by any finite number) once 
the origin is shifted. In case of the 80 -term IGRF 1965*0 model, Hilton 
and Schulz [ j shown that a 195-term representation (n 4 = 13 ) 

of the field in offset dipole coordinates maintains sufficient accuracy 
for practical purposes and have published its coefficients. It may 
also be noted here that coordinates of the type discussed here seem to 
be useful in describing Jupiter^s main magnetic field, as observed 
by the Pioneer 10 spacecraft |* Smith et al. , 197 ^^ • 


(f) PRACTICAL PROBLEMS RELATED TO UNIQUENESS 


The practical derivation of 7 for the main geomagnetic field 
involves many obseivational problems, such as the elimination from the 
data of effects of external sources and of magnetic anomalies in the 
earth's crust. Assuming these to be solved, the problem can be viewed 
as involving the derivation of 


7 = +- 7 ^ 


(21) 


where 7^ is some previously known model of the potential and 7^ is 
a relatively small correction which is to be derived from the given data. 
If 7^ is expanded as in (ll) in spherical harmonics involving some 
number N of unknown coefficients, then any observation related to the 
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field can fee reduced to a linear equation involving these coefficients. 
The relation may be exact, e.g* when ary component of B is observed, or 
it may be a linearized approximation , as happens when one observes the 
field's magnitude jBj (see later on) or any of the direction angles 
defining B » 

Thus the results of a world-wide magnetic survey tend to reduce to a 
large set of linear equations in K unknowns. The coefficients of 7 are 
then derived by the least squares procedure as those coming closest to 
fitting the set and, in case approximate linearization was used in obtai- 
ning the equations, 7^ is generally added to 7^ to give an improved 
initial approximation, after which the procedure is repeated one or more 
times [^ Cain et al. , 196 ? ] . 

An interesting problem related to such procedures is that of the 
uniqueness of the result r how can one be sure that a given data set 
leads to a unique choice of 7 ? 

One case in which this question has led to unexpected results 
involves the derivation of 7 from observations of jBj (commonly 
denoted in geomagnetism by F , a practice which will be followed here). 
Such observations are easily performed aboard spacecraft with sui accuracy 
of about one gamma j^ Cain , 1971 a^d references cited there a-nr^ they 
have the advantage of not requiring precise knowledge of the attitude 
of the sensor. If one neglects the altitude variation of the orbit, the 
uniqueness problem reduces to the question whether 7 is uniquely 
determined by observation of |7 7I = F over the surface of a sphere. 

Linearizing the expression of F obtained from (2I) gives 

and this can be used iteratively (as described before) to derive the 
coefficients of 7 , with 7 added to 7^ after each iteration step 

to provide a better starting approximation for the step following. In 
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computer simulations in which a known expansion (U) was used to generate 
F this method recovered the input coefficients quite accurately within 
steps. Consequently it was generally accepted and widely used in the 
treatment of actual data. 


Meanwhile^ however. Backus ^ 1968 , 1970 * 197^ investigated the problem 
mathematically o He first found that F and y were uniquely related when- 
ever the series (ll) was finite | Backus , 196 ^ , then he showed by actual 
counterexamples that this was no longer universally true if the series 
■ ( 11 ) was infinite f Backus , 1970~ | and finally he proved that this ambi- 
guity was removed if F was observed over a finite volume in space 
j ^Backus , 1974 j . 

At the same time evidence began to accumulate suggesting that models 
derived from F observed in near-circular orbits fit the vector field 
far less accurately than they fit the distribubion of the field's magnitude. 
In addition, different models derived from F exhibited relatively large 
differences and this prompted Hurwitz and Knapp ^1974] to conduct simulated 
recoveries similar to those described before but with data contaminated 
by finite "noise", as occurs in practical situations. In such cases y 
is not recovered exactly but finite errors remain and Hurwitz and Knapp 
found that the fit between the input vector field and that derived from 
the output model was decidedly inferior to the fit between input F and 
output F • 


Stern and Bredekamp ^ 197? J independently .obtained similar results and 
also showed that such effects were connected to the counterexamples of 
Backus . These counterexamples bear a special relation to the 

dipole field and the fact that the main geomagnetic field is dominated by 
its dipole component establishes a connection between them and the 
problem. In particular, it turns out that in the presence of finite 
"noise" certain sequences of harmonic terms can exhibit enhanced fluctu- 
ations which degrade the fit to vector data much more than they degrade 
the fit to F. Such enhanced fluctuations were in fact obtained in 
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computer simulations and also in analysis of some of the runs of Hurwitz 
Knapp I" 197^ ] and this strengthens the suspicion that some current 
models based on F are in fact less accurate than has been claimed. 

The uniqueness of y derived from scsne other types of data has 
been examined by Kono ^ 197 J in order to evaluate the correctness of 
some reconstructions of ancient geomagnetic fields. For instance^ he 
proved that magnetic declination observed at the earth’s surface does not 
in general define 7 uniquely, 

(3) TOROIDAJj A2® POLOIDAL VECTORS 

(a) INTRODUCTION 


A general vector field V may be represented by three scalars ^ 

and 4^ , in the form 

1 2 

I = (25J 

The advantage of this form is (as can be verified by carrying out the 
algebra) that if V satidfies the vector Helmholtz equation 

V +• V = 0 (24) 

then each of the 4^^ satisfies the appropriate scalar Helmholtz equation 

y. + k^ 4 ^. = 0 (25) 

(including the case k^= 0 , vhen (25) becomes Laplace’s equation). In a 
similar vayj if equation (24) Is written in cartesian coordinates, it is 
resolved into three scalar equations of form (25), one for each cartesian 
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component of V ; this is useful if the houndaiy conditions are easily 
expressed in cartesian coordinates, whereas (23) hears a similar relation 
to spherical coordinates. Only a few systems of coordinates allow such 
direct conversion of the vector equation to the corresponding scalar one 
C Senior , I960 j Morse and Feshhach , 1933 * chapt- 15 ] - 

Because (2^^.) arises naturally in wave propagation problems, the repre- 
sentation (25) first introduced in that context ^ Mie , I906 ; Hansen , 
1933 ] . If a magnetic field B is represented in this manner, is 

generally omitted, since it must then be a hannonic function and its 
contribution (as will be shown) may be represented by part of 4^^ o The 
representation using this form 

B ^ (26) 

was first introduced in connection with geomagnetic dynamo theory by 
Elsasser [ 19^3 , 19^6 , 19^7 * 1936 1 who called the two components the 
toroidal and pololdal components of B , respectively. 

These names are still used. To get some intuitive feeling for their 
significance, consider axisymmetrical fields, In which neither 4^^ 
nor 4^ depends on the azimuth angle 4* . A toroidal field of this 
type is then aligned with the 4 direction and has field lines 
circling the z axis, while the field lines of a poloidal axisymmetric field 
are orthogonal to the 4 direction and are confined to meridional planes, 
as another axisymmetrical example, consider a poloidal field with field- 
lines covering the surfaces of a family of nested toruses. When a toroidal 
component is added to this field, field lines acquire an azimuthal slant 
so that instead of staying in planes of constant 4 ‘they spiral around 
the axis of their torus, in a way resembling that found in toroidal 
plasma confinement devices. 

Some of the elementary applications of (26) to dynamo theory have 
been described by Elsasser . For instance, in a conducting 

sphere both components of (26) will tend to decay with time: differential 
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rotation - such as is observed on the sun and Jupiter, both of ‘which have 
magnetic fields - will amplify the toroidal component, but the strength 
of this process is proportional to the poloidal field which ultimately 
decays, unless it is continuously maintained by some feedback imechanism 
originating in the toroidal component, Parker showed, that radial 

flow in the rotating frame may induce cyclonic s'wirling which indeed 
leads to such feedback. Various theories of solar magnetism based on 
such considerations have been advanced £ Babcock , I96I ; Leighton , 196 9 J 
but the details are beyond the scope of this review. 


(b) VECTOR POTEKTIAL FOR A CTJRL-FREE FIELD 


Two identities useful in handling (26) are 

Vxy^£ = (27) 

V A V X 4 ^^ r = V - X C28) 

If B is curl-free and satisi'ies (2) , with y expanded as in (ll) 
and (20), then it can be represented as a poloidal field with ^ a 

2 

harmonic function simply related to 7 Stem , 196^ j . To see that, 
note that if 4^2 is a series of spherical harmonics, the same holds 
true for '^/Dr " Furthermore, the last term of (28) vanishes 

in that case, so that only the gradient of a harmonic function remains. 
This allows formal identification 

7 = - (ry ) (29) 

2 

and a term-by-terra comparison of expansion coefficients. Specifically, 
if 7 has expansion coefficients (g^, h^, h® ) as in (ll) and (20) 

and y has corresponding coefficients (cf , if , g" ) , then 
2 iT -o n’ n n n ^ 
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n 

11 

n 

= 




( 50 ) 

n 

= - 

n 

- - h^/(n^l) 
n 

To obtain the vector potential 

A for a 

given y . expanded as in 

ll) and ( 20 ), equations ( 30 ) are 

used to derive the expansion terms 

'f' j after which 
2 

one calculates 




A =7x4^ 
“ 2 

r 

(31) 


(c) MODELS OF THE MA.GKETOSPHERE 


Curl-free models of the far-away geomagnetic field, as were described 
in the section on the scalar potential, are incapable of great quantitative 
accuracy, because the region which they describe contains an appreciable 
current density. As an alternative. Mead and Fairfield 1975 j Fairfield 
and Mead , 1275 J tried to represent the field in such regions by 
expanding each component of B in powers of cartesian coordinates 


B. 


1 



n,r,s 


a. 

inrs 



n-r-s 

z 


(52) 


Potential fields such as those of (20) are easily reduced to this 
form but they comprise only a limited subclass of such fields. 


Two problems, arise, however: the first is how to assure the vanishing 
of ^ and the second is the difficulty in controlling the model 

near the boundary of the region for which it is derived, due to the pre 
of positive powers of x , y and z • 


To overcome the first difficulty Mead and Fairfield derived 
the relations between the coefficients of (^2) which are required in 
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order for y -B to vanish. They then analyzed magnetometer data from 
space and derived a "best fitting set" of coefficients for (32) by a 
least-squares procedure similar to the one used for deriving the main 
field, described in section (2-f) , but based on vector data. In that 
calculation the relations which assure the vanishing of y *B were 
introduced as constraints and were handled by the method of Lagrange's 
multipliers o 

A simpler method is however available: noting that 

and expanding ( 26 ) in cartesian coordinates shows that if 4^^ and 
4^ are expanded in the form (52), then the resulting components 
also have this form. In particular, if the largest power n of 
variables having the dimension of length is N for 4^^ and N-1 
for 4^^ , then Its value for the expansion of will be N . 

Kot all terms in the expansions of ^ produce independent 
contributions to B , because any part of which depends on . 

r alone does not affect -B . The terms with n - 0 belong to this 
class and should therefore be omitted. Consider next the three quad- 
ratic terms proportional to x^, and • These can be combined 
to give 5 independent terms proportional to X^, and and 

the last of these does not contribute to B, so that it can be omitted. 
Similar arguments show that for any integer k , one coefficient can 
be eliminated among those representing homogeneous polynomials of 
degree • 2k in the expansion of or y 
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One advantage of this approach is that no precautions are required 
to ensure the vanishing of 7 «>B • It is also easily generalized to 

overcome the second difficulty mentioned earlier, hy devising models 
in which ^ and are expanded in the foim 




a . 


ra S 

X y 


n-m-s 

z 


exp(-r/rQ) 




where is some chosen scale distance (in principle, several series 

with different choices of Tq could he used). It is more convenient 
to use in this case spherical harmonics 


V 

= z. 


n 


n,m, s 


exp(-r/r^) 


m 


[. 


Jnms 


cos s 


f + 


V 


jnms 


sin s'f 


( 35 ) 
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from -which the components of B are readily derived* It should he noted 
that the expansion of (34) is not equivalent to that of (3^)j iii fact, 
the latter expansion is equivalent to 




n 


r 



jnuvw 


u,v,w 


(x/r)^(y/i*)^(zA)^ exp(-r/r^) 


which includes the terms of (34) if ui-vi-w = n hut also contains 
additional expressions* Some comparisons between this method and that 
of Mead and Fairfield have been performed (D. Stern, unpublished) ; 
while the improvement is rather small compared to the inherent dispersion 
of the observational data, this approach should be useful in future 
studies by providing more flexible analytical expansions. 


Magnetospheric models somewhat similar to those described here have 
been used by Olson and Ffitzer Cisn] • However, the derivation of 
these models does not involve : it starts by fitting a system 

of currents to observed fields Olson , 1974 ] and then approximates 
the cartesian components of their fields by general expansions similar 
to ( 52 ) and ( 54)0 The resulting models are not automatically divergence- 
free • 

It should be stressed that there exists no assurance that "global" 
models, representing the entire magnetosphere by one series of analytic 
functions of a given form, can do so with arbitrary accuracy. The power 
series ( 32 ), for instance, is an ordinary Taylor expansion and can represent 
functions only in a certain neighborhood of the origin, where it 
converges: this region may not only be finite, it may be smaller than 
the magnetosphere, in which case any representation of the field in 
this way contains a finite error. The representation (54) appears at 
first sight to be more useful but lt,too, is subject to the same restrictions 
since it just represents the cartesian expansion of expCr/r^} . 
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At present^ observational inaccuracies are of the same order as those 
introduced by global models and such models are therefore in general 
use. Ultimately^ hoTwever, it may be necessary to adopt a different 
approach. Two extreme cases could then be considered - an extension 
of the global method by which the magnetosphere is divided into a number 
of large regions, each with its own expansion, and numerical represen- 
tation by means of a tabulated network of values (e.g. of It 

should be realized, however, that the numerical representation can also 
be viewed as a local expansion, by means of the interpolation formulas 
used, and by making the numerical grid more sparse and the interpolation 
formula more powerful, the gap between these two extremes could in 
principle be bridged. 


(d) SPHERICAL VECTOR HARMONICS 


The spherical harmonic functions 
- or, alternatively 


1^(0 ) cos ra^ 


and ( 9 ) sin m'f 
n 


n 



( 36 ) 


(in the notation of Morse and Feshbach [ 15^3 ] , p. I898-I9OO, which 
this section will adopt) - form a complete set in which arbitrary 
analytic functions of ( 0 ,^ ) can be expanded. Similarly spherical 
•vector harmonics are 3 sets of vector functions of 0 and 'f , denoted 
m ^ — n m ^ m ’ which find use in the expansion of vector 

functions. 


Many problems in wave propagation lead to a vector field V 
represented as in (23) and satisfying (sfe) for some values of k . 
This gives 3 independent equations of the form (25} and separation 
of the r variable shows then that ^ ^ have the form 




n.m 


A . F (r) ) 

nmi ran n * 


(37) 
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where F (r) stands for Hankel fimctions if k 0 and powers of r 
if k ~ 0 * ’When ( 5 ?) is then substituted in (25) it is found that all 
terms are expressed as the product of functions of r with some member 
of one of the following three families of spherical vector harmonics 


P 

— n^m 


—n,m 

C 

— n,m 


r x“(e,-f ) 

“• n ^ 

r [n(n + l)] Vx^Ce.'f) ( 38 ) 

[n(n-t-l)] V:<rx“(6,'f) 


Note that any dimensional dependence on r in the above definitions 
cancels out and also, by (27)^ that for any pair of values of the indices 
n and m the three above vectors ^re orthogonal. Other features of 
these vectors (and of their real and imaginary parts, since by (36) they 
represent complex quantities) are described by Morse and Feshbach [M2] • 

Many of the representation discussed earlier can be easily expressed 
in terms of spherical vector harmonics. For instance, the expansion ( 55 ) 
is formally the same as that of (57) and therefore the results of its 
substitution in (26) can be expressed in terms of the vectors listed in 
( 58 ). 
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(!].) EULER POTENTIALS 


(a) DEFINITION AND GENERAL PROPERTIES 


An intuitive description of a magnetic field, first introduced 'hy 
Michael Faraday, is provided hy magnetic field lines (or "lines of force" 
as he termed them) - lines which are everywhere tangential to B , Such a 
description is completely analogous to the representation of the velocity 
field of an incompressible fluid by means of streamlines.. 

For a long time field lines were mainly used as a device for visualizing 
the magnetic field. However, with increasing interest in particle motion 
and transport properties in a near-collisionless plasma - both in space and 
in the laboratory - field lines assumed new importance, since both these 
effects are channeled along them. This created interest in mathematical 
descriptions of the magnetic field which explicitely contain representation 
by field lines; such a description has been available in fluid dynamics 
and it involves two scalar functions ^ and ^ such that 

B = Voc X vp (39) 

The functions ) are generally known as Euler potentials 

because it was Leonhard Euler who first Introduced such a representation 
in the context of fluid dynamics; they appear to have been Introduced 
into plasma physics by Grad [ Northrop and Teller , I96O ; Gardner , I939 J , 
A short review of their history and properties was given by Stern [ 1970 ] 
and the reader is referred to that work for details extending beyond the 
present brief description. 

The basic properties of are easily derived from (39)* have 

= 0 

(40) 

B. Vp =0 
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This means that surfaces of constant ot and ^ are everywhere 
tangential to B - and this property extends to lines along which two 
of such surfaces intersect, which thus are field lines. Note that (59^ 
implies (40), hut not vice versa. Two functions (u, v) with the property 
of (40) 

B»Vn — B®Vv ” 0 

do not satisfy (39^ j "the general case, but rather 


B = w ( 7u X Vv ) 

where w is an arbitrary function of (u, v). Functions such as (u, v) 
may be called unmatched Euler potentials and have been introduced by 
Sweet [ 1930 ] and Dungey [ 195^ 1 • They are conserved along field 

lines and are therefore functions of ( 0 ^, p) • 

In general (subject to restrictions of uniqueness and single valuedness 
noted later) field lines of a given configuration form a two-parameter 
family, in which the field line corresponding to the values ( 
of the parameters are represented by the intersection of the surfaces 


o( (x, y, 2 ) = 

^ (x, y, 2 ) = fo 


(if-l) 


As an example, a simple model of the geomagnetic field is given by 
an axial dipole of moment ; if (as in eq„ ll) a represents the 

earth’s radius, one convenient choice of ^ in spherical 

coordinates (r, 0 , , is given by 


oC = a (a/r) sin^ 0 (42-a) 

P = at 


(42-b) ■ 
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Given any two numerical values { Po> of ^ and j a field line 
of this field is defined by (^l). Note that (i|2-b) indicates that field 
lines lie in constant meridional planes, while (42-a) gives the shape 
of a field line within such a plane. Other examples will be presented later, 

Euler potentials are not uniquely defined: for Instance, ^ can be 
incremented by any function of p without ( 59 ) being violated. More 
generally, C f ) may be replaced by (o<', p’) , provided the new 

Euler potentials are functions of the old ones and the Jacobian of the 
transformation is unity: 

p')/^(o<,p) = 1 (i,5) 

Given a magnetic field represented as in (40), it is generally possible 
to choose for ' (at least inside a restricted region in space) any 
’w®ll“behaved function and then derive a corresponding 

conjugate Euler potential p ' • An application of this property will be 
presented at the end of this section, 

A vector potential corresponding to ( 59 ) is 

A = o<vp (44) 


and it has the property A»B = 0 , 

With all these useful properties, there exist two important drawbacks 
to Euler potentials which limit their application. The first is non-linearity: 
the representation ( 59 ) is non-linear, since it involves products of the 
derivatives of oC and ^ • This means that superposition does not hold 
for Euler potentials: given the Euler potentials for two fields and 

^ , the functions for the combined field does not equal the 

sum of the Euler potentials expressing and separately and, indeed, 

might be difficult to derive, even when the Euler potentials of the com- 
ponent fields are known. In practice this greatly limits the class of 
fields for which analytic forms or even analytic approximations of 
are readily derived. 
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The other dra'w'back, more Importaxit in the laboratory than in space 
applications,, is the possibility that the labeling of field lines by 
C ot, p ) is not single valuedo In toroidal confinement devices, for 
instance (tokamaks, stellarators), if a field line is labeled by ) 

in some limited region of the field and is then followed outside this 
region, it may ultimately’’ enter the region once more. In fact, such a 
field line usually returns to the same region again and again and with 
each return it coincides with some previously labeled field line - in 
general, one with different values of (o<, In such cases the labeling 
is single valued only if a. limited region is considered and is not possible 
when the entire field is represented. For further details about these and 
other properties of Euler potentials the reader is referred to the review 
Stern I97O "j , 

As an illustration of an application consider the motion of trapped 
particles in the magnetosphere under the combined influence of a static 
electric field 

E = - 

and the geomagnetic field B . Because of the high conductivity along 
field lines the electric field (except in some special regions, perhaps) 
comes close to being orthogonal to B , i,e. 

= 0 (lf5) 

By (4o) this condition is neatly expressed by the requirement that 0 
depends only on ( , P ) 

( 46 ) 

A charged particle of very low energy, conserving the two lowest 
adiabatic invariants but with negligible magnetic drift, will move in 
this field with the drift velocity 

Vg = (S»V 4 ’)/b^ (47) 



- 50 - 


From this v <> V</> = 0 and therefore the particle stays on a single 

“E 

equipotential surface throughout its motion. 

Now as stated earlier, any "-well-hehaved” function oi ' of 
can be chosen to replace oi and an appropriate "conjugate potential" ^ ' 
can then be found for it. In particular, let the role of oi ’ be played 
by and let the conjugate potential be denoted 

i ,e . 

B = 7^ X (48) 


It is instructive to derive the rate at which 'V changes at the 
location of the drifting particle. By (4-7) and (48) 

d^f/dt = v_.Vy - 1 (49-a) 

““Xi 

i.e , 

^ ■ ■’^0 (49-b) 

Thus a swarm of particles starting from a surface of constant at 
t = 0 will always share the same value of , The averaged motion of 
such particles - that is, the gradual change of their guiding field lines - 
is conveniently studied in the ) plane, where to every pair of 

values of the Euler potentials - and hence, to every possible guiding 
field line - there corresponds one point (in the geomagnetic field the 
equatorial plane can play such a role)«In this plane the lines 

^ = constant 

can be viewed as rays along which very low energy particles propagate, 
while lines of constant 4^ act as "wavefronts" marking the particles' 
progress. Figure (l), taken from Stern J gives lines of constant 

^ (solid) and 4^ (broken) for some particular geomagnetic model, in 
which the earth’s field is approximated by a dipole field and some analy- 
tical model of ^ is assumed. 



- 51 - 


The preceding is readily extended to paarticles of finite energy 
conserving the magnetic moment M and the second invariant J , 
except that now ^ is replaced hy the averaged guiding center 
hamiltonian K( ^ j) introduced by Northrop and Teller I96Q 
This generalization has been developed for the dipole field by Chen 
and Stern ^ 19'^ j 3 "who provided an analytical approximation for K 
and expressed its conjugate potential by means of numerical Integration. 


(b) EXAMPLES 


In general, simple analytical forms for the Euler potentials can only 
be found for classes of fields vith some type of symmetry. For instance, 
two-dimensional fields of the form 

B = B^(x, y) X + B^(x, y) J (50) 

have Euler potentials 

o< = / dy + f (x) = - / dx t g(y) 

(51) 

^ = z 

where the functions f(x) and g(y) are obtained by deriving ©4 in 
two different ways (because V = 0 , two such equivalent derivations 
are possible). The hamiltonian for a particle with mass m and charge q 
moving in such a field can then be written 

H = (l/2m) + C - q P<(x, y) (59) 

Since p is a constant of the mction it follows that the motion 
z 

reduces to that of a two-dimensional free particle in a non-negative 
potential 


V 


^ Pz - 


q 0< )^/ 2m 


(55) 
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For all such motions, lines of constant oC in the (x, y) plane - 
which trace the pattern of field lines ~ are equipotentials of V ® 

The value of V attached to ai^ such line and the general configuration 
of regions of high and lov potential depends on the constant p appro- 
priate for the given motion« One interesting configuration of this class 
occurs in the neighborhood of an X-type neutral line ^ Russbridge , 1971 ^ 
other examples and generalizations have been examined by Stern ■ 

Axisymmetricai poloidal fields, by (26), can be represented by Euler 
potentials with ^ = 'f , since 


V ^ V X r f.(r, e ) 


V X 'p ( /-ae ) 


= - v[r sme ^ vy 


( 54 ) 


The dipole example of eqs. (ii-2) belongs to this class. It is also 
possible to derive simple Euler potentials for pure toroidal fields: 


V X r - VH'i X Vr (55) 

However (as noted earlier) when fields of both kinds are combined 
no simple way of deriving Euler potentials exists. 

A case of practical interest involves the main geomagnetic field, 
usually represented by a scalar potential 7 as in eq, (ll). No 
ana3^ical method is known by which (o(, may be simply derived 
for such a field, but a perturbation technique is available for 
obtaining them approximately. 

The method is based on the observed property that the dipole 
component in (ll) dominates the expansion, especially if 7 

is expressed in dipole cuordlnates » In such coordinates let 
be the dipole Euler potentials of eqs« (te) (the subscripts zero refer 
here to zero-order approximations and are unrelated to those appearing 
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in equations ), 'while 7^ is the corresponding scalar potential. Let 
7 ) he the first-order corrections, linear in the higher 
harmonic coefficients f applied to these three quantities. 

Substitution in (39) then gives 




from which 


Vo(^-V7o = 
Vf^‘V7Q = 


(57) 


Expressing { , fl , 7 ) in terms of ( 0 0 ) allows these 

1 1 O' r 0' 

equations to be integrated and provides ( , 0 ) within arbitrary 
functions of (o^^, p ^). However, (57} contains onl^y two of the three 
scalar relations implied by the vector equation (56); if the remaining 
scalar relation is now invoked, the arbitrary added function may be 
determined. Details have been derived by Stern [mi] and related 
expressions were first obtained in a different context by Pennington 
|^ 1961 » 1967 ^ • The calculation has not been extended to higher 

order, except for one effort by Hassit ^ 196 in which second-order 
unmatched Euler potentials were derived. 


The same perturbation method can also be applied to the external 

harmonics of section (2-i) and in particular to the simple model 

described there, with the expansion of 7 involving only the three 
0 _0 

coefficients g^, and gr . As might be expected, the resulting 
expansion breaks down at large distances, but it also becomes unrealis- 
tic near the z axis. It provides a fairly good representation of the 
field in the regions occupied by trapped particles; in the noon- 
midnight cross section of the ) model^n the region in which the 

model field departs markedly from the curl-free field from which it 
was derived, there exists the added feature (not found in the curl-free 
model) of a neutral point on the night side of the dipole, in addition 
to two such points on the day side, corresponding to the polar cusps. 
Further details have been given by Stem I967 J and Kosik ^ 1971 a , 
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A model vith some similar properties (tut much smaller deformation) has 
been developed by Parker [I960] - 

While accurate representations of the magnetospheric field by means of 
are difficult to obtain and no iterative procedures for deriving 
them have been published, it is relatively easy to devise simple models 
vhich exhibit appropriate qualitative properties. For instance, a distor- 
ted dipole field, compressed on one side and stretched out on the other, 
is given by 






sin^0 / r ^ k sin®G (1 - cos'f’ 

a? 


(58) 


where k is an adjustable constant. It is easy to distort this model 
further, in such a way that all its field lines are contained -within 
the surface 

f(r,e,'f) = 1 (59) 

This is accomplished if o( is multiplied hy ( 1 - f ) ; in order that 
the field will continue to approximate the dipole near the origin, f should 
tend to zero as r 0 . For instance, if the field is to be contained 

inside the paraboloid 


r = ^q /(1 + sin 0 cos'"P) 
one can multiply o( of (58) by the factor 
1 


(r/r^) (1 -t sln& cos'f ) 


m 


(60-1 


( 60 -b) 


Increasing the adjustable parameter m narrows down the transition 
zone between the main field and the boundary field. Further improvement can 
be obtained by shifting the focus of the paraboloid away from the origin, 
but it should be noted that in all these models the boun(Jary is always 
connected to the two dipole field lines which emerge along the z axis. 


As a final note, axisymmetric models of this sort have been construc- 
ted by Parish et al . £ 1974 j for the magnetic field of Jupiter, following 
its observation by Pioneer 10 Salth et al. , IST^ J . 
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(c) THE L PAHAMETEB 


Before the ) system for the magnetosphere was developed an 

alternative approach to the labeling of magnetic field lines was 
introduced in the form of the so-callcd L parameter due to Mcllwain 

[ 1961 . 1966 ] . 

As an intuitively meaningful quantity the L parameter is e^ctremely 
useful: it has dimensions of length in units of earth radii and in a 

S 

dipole field its value at any given point equals the maximum distance 
from the dipole attained by the field line passing the given point. In 
a perturbed dipole field (and the parameter is only defined for dipole 
fields and perturbed dipole fields) the value b(r) at some given 
position r still approximately equals the largest distance from the 
origin which is attained by the field line through r • However, the 
generalization of L to perturbed dipole fields is not performed in a 
way which aims at presei*ving the equatorial crossing distance. Instead, 
the definition tries to preserve a different property, useful in studies 
of trapped radiation - that if a charged particle is trapped along any 
field line with a given value of L, as it gradually drifts around the 
dipole then all field lines which it passes will have the same value of 
L . 


This property is only approximately achieved. Indeed, it cannot in 
general hold for all particle orbits In an asymmetric field, since two 
particles starting on a given field line in such a field with different 
mirroring points will in general follow slightly different surfaces ("drift shells" 
during their drift (this is known as shell splitting). However, for the 
geomagnetic field the separation between such surfaces tends to be small 
and consequently the values of L encountered by a particle in its drift 
generally vary by less than 2 % . 

In this section the relation between L and ) will be traced, 

the inherent limitations of L will be clarified and suggestions will be 
described for ways in which the concept can be generalized. 
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In a static near-dipole magnetic field, trapped particles tend to 
preserve the tvo lovest adiahatic invariants (the third invariant is not 
needed here ) 


M- = /2niB 


J = ^ P, dX (6l) 

h 


where m is the rest mass, (p^ , ) are the components of the 

momentum parallel and perpendicular to B and integration is along 
a field line. Since the magnitude p of the momentum is also conserved 
(in the absence of electric fields, which will, now he assumed) it is possible 
in this case to replace (n, J) by two related quantities I) which 

are also constants of the motion but do not depend on the particle’s 
energy: 




I = J/£P 


/ ® (1 - b/b dl. 
B' “ 


( 62 ) 


where the integration of I extends between mirror points at which 
B = » As the particle drifts from one guiding field line to the next, 

it always chooses that one of the adjacent field lines on which the 
value of I , evaluated between fixed values of B , is the same. 


As the particle drifts its guiding field lines gradually trace a 
surface - a closed surface for trapped particles - called a drift shell 
or a drift surface. Because drift shells are tangential to B their 
equations have the form 

= constant (65) 

For each pair of parameters (l, there exists a drift shell and 

therefore the entire collection of drift shells can be represented by 
some function F as 
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= constant (64) 

(the constant may be chosen as 0 or 1 vithout loss of generality) » 
Isolating gives 

ck = G(l, p ) (65) 

In a dipole field ^ may be eliminated because of axial symmetry; in 
that case, if the Euler potentials defined in (^2), then 

(65) reduces to 

*^0 = 

Kov (l, B^) can be replaced as parameters characterizing drift shells 
by any tvo "well-behaved functions depending on them j^ Stern , 1968 *] • In 
particular, in the case of the perturbed dipole field, 3-^7 

function L(Gq) can be used as one such parameter, while B^ may be retained 
as the second one. This characterization will have the advantage that in 
the dipole limit - due to (66) - the shell equation depends only on a single 
parameter or, L(Gq)_, not on two» In that limit all the trajectories 

starting from the saiiie initial field line trace the same surface and have 
identical values of G^ or L(Gq); by way of contrast, if such shells 
were classified by using (l, B^) , a finite range in both these parameters 
would be required. 

In a perturbed dipole field the shell equation is given by (66) 
and no single-parameter description is possible. However, if the pertur- 
bation is small, equation (65) will differ from (66) only by small 
correction terms. Denoting such terms by subscript "l” and using (Gq, B^) 
as parameters brings (65) to the form 

oi 




(67) 
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Thus the shell's equation depends strongly on and only veakly, 

through the correction term G , on the second parameter B » If G, 

1 Dl 1 

is ignored - or "better, replaced hy <«x> , its average over - 

then all particles characterized hy some value of G^ will stay close 
to the surface 


Cy = + <G^(G„,p> 


( 68 ) 


and thus their drift shell will still he approximately characterized by 


the single parameter • All the preceding also holds if 


G 


0 


replaced by L(Gq); in particular, Mcllwain's definition of L ^ Mcllwain , 
120] is equivalent to 


L(I, B ) 
■’ m 


^ b_^) 


(69) 


In Mcllwain's work the above function is approximated in two steps, 
beginning with the definition of auxiliary functions 


Y = In (L?Byg° - 1 ) 

X = In^ where ^ 


(70) 


The relation between X and Y is then approximated by a polynomial 

N 


n=0 

In the original derivation Mcllwain , 1961^ sixth-order polynomials 
were used and different expansions were chosen for each of 5 ranges 
of X . Later on Mcllwain , I966 J an improved approximation was derived, 
with K increased to 9 the number of ranges to 6 , and an inverse 

expansion, expressing X in terms of Y , was also provided. A much 
simpler and fairly accurate approximation has been proposed by Hilt on 
[1211] , who matched the analytical behavior of L at the lim its 


1=0 and B 


m 


oO with that of the approximation. His result is 
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lX/s° = 1 + + =^3^ (72) 

vhe re 

a = 3 (2)^/^/n « 1.5507 

1 

= [a + (3)^/'^ln(2 -I- 3^/®)] 

0 . 0475^^55 

a 3 - 0.456576 

The relative inaccuracy dL/L of this approximation is 10 "“^ or less. 


The L parameter is videly used in labeling field lines, in the 
folloving manner. Given a point P = (x, y, z) in the earth's field 
(■which is assumed to be given by the expansion (ll) of the scalar 
potential) , the integral I is derived for particles mirroring at P , 
by numerical integration. Using one of the formulas described earlier, 

L(I, B^) is derived for these particles and its value is regarded as 
a labeling parameter for the field line passing P • Of course, what has 
been derived here is an approximation to a g^/oCCx, y, z) where 0<1 is 
the perturbed-dipole generalization of <Xq in ( 4 £-a) (in principle a S^/o< 
could be chosen as a new Euler potential c< ' , but the form of ^ ’ is 

then complicated). The derivation of this approximation is encumbered by 
the need for a numerical integration expressing I , although this 
integration has been simplified - as far as computer use is 

concerned - by G. Kluge (ESRO Internal' Note 66 , October I97O) ^ho first 
transformed it by inversion with respect to the origin. 


This correspondence between L and works both ways and Stern 

^ 1968 j has used approximations of oL for perturbed dipole fields to 
obtain approximations to L(I, avoiding the need for deriving I 

by line integration. This method is already implicit in the work of 
Pennington I96I , I967 who derived the equations of drift shells in 


a perturbed dipole field by a strai^tforward perturbation method. 



- 40 - 


poloidal 

The L parameter can be generalized to any near-axisyrametric ^magnetic 
field : if (65) represents the family of drift shells in such a field, it 
is only necessary to derive the corresponding equation ( 66 ) for the 
underlying axisymmetrlc field and then choose or some function 

of Gq as the L parameter*. 


Such an approach was applied by Stem [m§] to the 5 -parametei' model 
of the magnetosphere discussed in the preceding section. At moderate dis- 
tances from the origin this field may be viewed as a perturbation of the 
axisymmetrical field given by 



(73) 


By a perturbation method Stern derived an approximate form for the 
equation corresponding to ( 66 ) in this field 


= cyi, Bj 

and denoted L' = a S^/g' as the generalized L- parameter. In principle 
the same generalization should be feasible for talcing into account the 
axisymmetrical ring current of the earth, for treatment of particle motion 
in slightly asymmetrical mirror machines and for similar applications. 


A number of quantities related to L are in general use. Among the 
best known of these is the invariant latitude f\ : given a point with 
some value of L , its "invariant latitude" satisfies 

cos^A = i/l (t 4 ) 

The underlying idea is that if all points in the perturbed field were 
to be mapped to a dipole field in a way that L (or o( ) were conserved, 
then A would be the latitude at which the field line through the point 
met the earth's surface. A similar generalization for the longitude is 
provided by the magnetic longitude 


Mcllwain , I966 ^ 



(5 ) LOCAL REPRESEKTATIOWS 


It is sometimes required to represent B in the immediate neighbor- 
hood of some given point P vithout regard to the field's configuration 
in the rest of space » Two important cases in which this happens are in 
the guiding center motion of charged particles and in plasma dynamics 
near neutral points: it is useful in such cases to introduce so-called 
local representations of B which are valid only in a restricted 
neighborhood of P » 

The most common local representation and the only one discussed 
here - is the Taylor expansion 

B(r) = Bq -h £.VBq + |rr;v7BQ+ ... (75) 

where all quantities with subscript zero are evaluated at the point P , 
regarded as the origin for r » Such a representation is fully equivalent 
to ( 52 ) and is best handled by the methods developed in section (5"c)^ 
i.e. by expressing B in the form 


B = 

V X 

r + V ^ V 

(76) 


and expanding 4^ and 

2 

in polynomials in (x, y, z). Let 


denote the sum of those 

terms 

in the expansion 

of y ^ which form 

homogeneous polynomials 

of degree J in (x, y, 

, z). Then it is 

easily 

seen that 





ltd 

c 

It 


V A V X r 



= 


(2) 

V X V A 


(77) 

~TT :7VBq = 


(3) 

V X 

^(2) 

4* Vxry^ 



and so forth, with the degree of the terms on both sides increasing 



^ h-2 - 


hy unity with each succeeding line. The coefficients on either side are not 
in general independent. On the left side the terms are related ly virtue 
of the condition 7 = 0 ; for instance, V is required to have 

zero trace, so that on3y 8 of its coefficients can he independently 
specified. On the right any combination of terras which depends on r alone 
does not contribute to the field and should be eliminated by the methods 
of section (5-c). The vector potential follows at once from (7?) j it should 
be noted that no simple relation exists between the expansion (75) of B 
and a similar one for A . 


In principle ^ Morse and Feshbach , I955J an alternative local 
expansion is possible in which r in ( 76 ) is replaced by the constant 


unit vector 


A 

z 


B = (T8) 

- -^1 “ *^2 ^ 


This, however, introduces a considerable ambiguity. In ( 76 ) no compo- 
nent of B can be both poloidal and toroidal, so that the assignment of 

such components to and ^ is never arbitrary. The analogous 

1 2 

statement does not hold for ( 78 ) and because of this, except perhaps for 
special cases, this expansion is not recommended. 


Local Euler potentials are similarly expressed as polynomials of ascen- 
ding orders 


with 


= 

^ 4 - 

o< + 


( 2 ) 


+ 


(79) 

(2) 

2o = 



£-VBo = 

Vo< * V| 

^( 2 ) 



(80) 


and so forth. The calculation here is greatly simplified if the coordinate 
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axes are suitalDly chosen, vith the z axis along and 


ot = ^ 


= V 


If then 

(2) 




= aiiX^ + + asa^^ + 2b22)xy + a^3xz + 


f - h^^x^ + + ^33^^ +. (h^- + higxz + h^3yz 


(81) 


it is found that the 4 coefficients a , h,. and h affect 

only orders higher than the first in ( 77 ) > v>hile the remaining 8 
coefficients are uniquely related to the 8 independent terras which 
determine V Bq « 

An interesting local representation, equivalent to (75) ^iP io the 
order of > has "been described by Siainbis and Korthrop [1266] . 

.A A A 

At a point P in the field let the unit vectors vL, M, NJ form a 

A 

right-handed orthogonal system with the origin at P , L parallel 
to B and M pointing towards the center of curvature of the field line 
through P « In this system the components of VBq can be expressed 
in terms of curvature and shear coefficients (loco cit^, appendix) 
while the coefficients themselves can be expressed in terms of the spatial 
variation of the basic unit vectors following a shift in P . Some 
applications of such a system of local vectors exist in plasma theory, 
e«g„ in the work of Wilson [1212] . 


The Taylor expansion i'^) forms a focal point of the theory of neutral 
points I^ Dungey , Igs] - A first-order neutral point exists at F 

of vanishes but V does not, so that to lowest order, in the 


vicinity of P 
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B (82) 

Even though B vanishes at P there may exist singular field lines 
which pass through P and at sufficiently small distances from P such 
field lines can then be approximated by straight lines. In the region of 
interest the magnetic field B depends linearly on (x_, y^ z) as 
measured from P and therefore, if | is the radial distance measured 
from P along one of the singular lines, its value on that line should 
satisfy (to lowest order) 

B = (85) 

where is some constant. By (^) 

= I'VBq (84) 

and therefore X must be a real eigenvalue of . There may exist 

either 1 or 5 such eigenvalues ; accordingly, there will exist 1 
or 5 such singular field lines and P will be classified as an 0-type 
or X-type neutral point, so named because field line configurations 
near P resemble either the letter 0 or the letter X (Figure 2) • 

A great amount of theoretical work and of interest has focused on 
X“type neutral points, since they play a central role in processes 
which change the topology of magnetic field lines; for a review the 
reader is referred to the work of Vasyliunas 197^ j • Surprisingly, 

only limited experimental work exists on the subject ['e.g. Baum et al. , 
2m] > nor has there been much interest in magnetic field configurations 


in space with 0-type neutral points 
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APIERWQRD 


This is a draft of a review article and additional relevant material 
may be added to the final versiono Hie author welcomes any suggestions 
for such additions as well as comments of a general nature o 

Two additions which will be incorporated concern the description 
of two extraterrestrial magnetic fields by harmonic expansions similar 
to eqo (ll). Altschuler and Kewhlrk. j] 1969 ] applied such an expansion 
to the magnetic field of the sun and their work was extended by Schatten 
[ ^971 ; Schatten and Howell , 1971 1 • In addition, Acuna and Hess 
[ 1^75 ] extracted a model of the external field of Jupiter from the 
observations of Pioneer 11, which passed close enough to the planet to 
allow the larger non-dipole harmonic terras to be estimated. 

Another such addition concerns the work of Voigt ^ 1972 ^ j who 

developed a general class of current-free magnetospheric models. In 
Voigt’s approach the bounding surface (i.e, magnetopause) is deduced 
from observations and serves as a boundary condition to which an 
expansion of 7 is then fitted. 
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CAHEIOWS 


Figure 1 


Figures 2 


TO FIGURES 


- Lines of constant electrical potential <f> { Oi., (solid) 
and of constant conjugate potential ) (dashed) for 

a dipole magnetic field in the equatorial plane, using a 
simple analytical model of the earth's electric field* 


- The behavior of magnetic field lines near an X-type (2-a) 
and an 0-type (2-b) neutral point* In Figure (2-b) the 
pattern belov the neutral point mirrors the one above it 
and the main axes need not be orthogonal* 
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